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Our knowledge on the nucleon spin structure has greatly improved over the last 
twenty years or so, but still many fundamental questions remain unsolved. I will 
try to review some of the puzzling aspects of the origin of the nucleon spin. I 
will emphasize the connection with several sum rules and, when using this tool, 
the relevance of some kinematic regions for testing them in the QCD dynamics 
framework. 



1. Introduction 

The reply to the question " why do we need spin in high energy particle 
physics? ", on very general grounds, is twofold: first we want to learn 
about hadron structure and second we want to test perturbative QCD in 
the spin sector. More specifically concerning the nucleon structure we need 
to determine the unpolarized parton distributions fN(x,Q 2 ) (N = p, n 
for proton and neutron), where / stands for quarks (u,d,s, ..), antiquaries 
(u, d, s, ..) or gluons. We also need to know the corresponding helicity dis- 
tributions Afjy(x,Q 2 ) and transversity distributions. All these distribu- 
tions are functions of the scaling variable x and the Q 2 dependence, so the 
scaling violations predicted by QCD, must be compared to experimental 
results. The data allow to extract various structure functions, unpolar- 
ized -Ff' n , -F3, .. and polarized g\' , from deep inelastic scattering (D1S), 
which are expressed in terms of the parton distributions. The relevant ex- 
periments are currently performed at CERN (Compass), DESY (Hermes), 
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JcfTLab and SLAC, but one should also mention the new RHIC spin pro- 
gramme at the polarized pp collider at BNL 1 , which will go beyond testing 
the QCD scaling violations of the polarized parton distributions. 

2. Many sum rules on structure functions 

There exists a number of sum rules for unpolarized and polarized structure 
functions, some of which are rigorous results and other which rely on more 
or less well justified assumptions. Let us first consider the charged current 
structure functions in neutrino DIS, for which we have two rigourous results, 
namely the Adler sum rule 2 (ASR) 

jf 1 g - F?{x)} = N U -N d = l (1) 

because 

N u = dxu va i(x) = 2 and N d = / dx d val (x) = 1 , (2) 
Jo Jo 

and the Gross-Llewellyn Smith sum rule 3 (GLSSR) 
f 1 dr 

Igls= ^[xF» p (x)+xF» p (x)]=N u + N d = 3. (3) 
Jo ^ x 

The ASR is exact and receives no QCD corrections, but its experimental 
verification is at a very low level of accuracy 4 . The GLSSR gets a negative 
QCD correction and the CCFR data 5 gives I GLS = 2.55 ± 0.06 ± 0.1 at 
Q 2 = 3GeV 2 , in fair agreement with the theoretical prediction. 

Next, let us consider the unpolarized electromagnetic structure func- 
tions for proton and neutron F% ,n (x) , for which we have the Gottfried sum 
rule 6 (GSR). If one assumes an SU(2) symmetric sea, i.e. u(x) — d(x), one 
can easily show, using Eq.(2), that 

Ig= [ 1 ^[FZ(x)-F?(x)]= 1 -. (4) 
Jo x J 

In fact the NMC experiment 7 has observed a large defect of the GSR, 
since their measurement gives at Q 2 = 4GeV 2 , Iq = 0.235 ± 0.026. This 
flavor symmetry breaking, more precisely d > u, is a consequence of the 
Pauli exclusion principle which favors dd pairs with respect to uu pairs, 
since the proton contains two u quarks and only one d quark. 

If we turn to polarized structure functions, there is first a fundamental 
result called the Bjorken sum rule 8 (BSR). It was derived about thirty 
years ago in the framework of quark current algebra and it relates the first 
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moment of the difference between g\{x) for the proton and g™ (x) for the 
neutron and the neutron /3-decay axial coupling 

£dx\£(x)-g?(x)] = ±g A /g v , (5) 

where gA/gv = 1-2573 ± 0.0028 is very accurately known. The BSR gets 
also a negative QCD correction and we will come back later to the test of 
this firm prediction of QCD. One can also derive sum rules for g\ and g" 
separately. These are the Ellis-Jaffe sum rules 9 (EJSR) which read 

r 1 i 

r?= / dxg%(x) = —(9F-D + 6As) and 
Jo 18 

f 1 1 

r? = jf dxg?(x) = —(6F-4D + 6 As) , (6) 

where F = 0.459 ± 0.008 and D = 0.798 ± 0.008 are the /3-decay axial cou- 
pling constants of the baryon octet and As = J Q As(x)dx is the total polar- 
ization of the proton carried by the strange quarks. One recovers Eq. (5) by 
taking the difference because F + D = gA/gv- 111 their original work, Ellis 
and Jaffc made the critical assumption As = 0, which allows to make defi- 
nite predictions for and T™. One gets = 0.19, in strong disagreement 
with the sixteen years old EMC data 10 F{ = 0.112±0.009±0.19. This large 
defect, which still remains, was first attributed to a large As, but this naive 
interpretation has been ruled out since then, on experimental grounds. One 
should keep in mind that to test all the above sum rules, one requieres an 
accurate determination of the corresponding structure functions in the full 
kinematic range < x < 1, which is never achieved in current experiments. 
Significant progress will have to be achieved, to reduce the missing region, 
where some dangerous extrapolations are needed. 

Concerning the other polarized structure function g2(x), which is related 
to transverse polarization, but has no simple interpretation in the parton 
model, it is possible to derive a superconvergence relation by considering 
the asymptotic behavior of a particular virtual Compton hclicity amplitude. 
This leads to the Burkhardt-Cottingham sum rule 11 (BCSR) 

/ dxg p 2 (x) = and f dxg%(x)=0, (7) 
Jo Jo 

for proton and neutron and from this result, it has been naively argued that 

gi{x) vanishes identically. In fact, alternatively a simple relation between 

gi and gi can be expected, namely g\(x) + giix) = 0. However 52(2;) is 

more complicated than that 12 and only part of it (its twist- 2 contribution) 
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is entirely related to g\ (x) , by means of the Wandzura- Wilczek sum rule 13 
(WWSR) which reads for J > 1 

0. (8) 



f dxx- 1 - 1 
Jo 



1 1 Mx)+ 9 r w (x) 



j 

Clearly for J = 1 one recovers the BCSR Eq. (7) and for J = 2 one has 

9^ W {x) = -9i{x)+ fgiiv)^- (9) 

J X V 

Finally, using the spin-dependent photoabsorption cross sections 
<T i/2(3/2) { v )i Gerasimov-Drell-Hearn 14 (GDH) have derived the following 
celebrated sum rule valid for real photons 

f°° r /x , ^dv 2ir 2 a 2 ,„„x 

Igdh= Wi/2(v) - tT 3 /2(v)]— = — rp-K . (10) 

where v is the photon energy in the target rest frame, v t hr is the pion 
production threshold, k the anomalous nucleon magnetic moment, M the 
nucleon mass and a the fine structure coupling constant. However the GDH 
integral can be generalized to the case of absorption of polarized transverse 
virtual photons with Q 2 

Igdh(Q 2 )= ^ [<t 1/2 (^Q 2 )-t 3/2 (u,Q 2 )}-. (11) 

One can show that, within a good approximation, one has in the scaling 
limit 

Igdh(Q 2 ) = ^T 1 , (12) 

where ri = J dxgi(x), so the GDH sum rule is connected to polarized DIS. 
The <5 2 -dependencc of the generalized GDH sum rule has been measured 
very accurately recently, with different targets, and we will come back in 
more details to some theoretical understanding of these data. 



3. The statistical approach for polarized parton 
distributions 

DIS of leptons on hadrons has been extensively studied, over the last twenty 
years or so, both theoretically and experimentally, to extract the polarized 
parton distributions of the nucleon. As it is well known, the unpolarizcd 
light quarks (u, d) distributions are fairly well determined. Moreover, the 
data exhibit a clear evidence for a flavor-asymmetric light sea, i.e. d > u 
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as mentioned above, and large uncertainties still persist for the gluon (G) 
and the heavy quarks (s, c) distributions. The corresponding polarized 
gluon and s quark distributions (AG, As) are badly constrained and we 
just begin to uncover a flavor asymmetry, for the corresponding polarized 
light sea, namely Au ^ Ad. Whereas the signs of the polarized light quarks 
distributions are essentially well established, Au > and Ad < 0, this is 
not the case for Au and Ad. Here we briefly recall how we construct a 
complete set of polarized parton (all flavor quarks, antiquarks and gluon) 
distributions. Our motivation is to use the statistical approach 15,16 to 
build up : qi, Aq^, q~i, Aq~i, G and AG, in terms of a very small number 
of free parameters. A flavor separation for the unpolarized and polarized 
light sea is automatically achieved in a way dictated by our approach. 

The existence of the correlation, broader shape higher first moment, 
suggested by the Pauli principle, has inspired the introduction of Fermi- 
Dirac (Bose-Einstein) functions for the quark (gluon) distributions 17 . After 
many years of research, we recently proposed 15 , at the input scale Qq = 
4GeV 2 



xu+(x, Ql) - _ ^ ^ + i + cMx/ _ ) + 1 , (13) 

TV~(r (l 2 ) - ^( X 0») lx2b I Axb ( U ) 

xu (x, y ) cxp[(a; + x + ^ + i + cMx/ _ ) + 1 , U4j 

Ar<X~ b+1 

xG( Xl Ql) = G , (15) 

exp(x/x) — 1 

and similar expressions for the other light quarks (u~,d + and d~) and 
their antiparticles. We assumed AG(x, Q§) = and the strange quark 
distributions s(x,Qq) and As(x,Qq) arc simply related 15 to q(x, Qq) and 
Aq(x, Qq), for q = u,d. A peculiar aspect of this approach, is that it 
solves the problem of desentangling the q and q contribution through the 
relationship 15 

*o + u + *olz = 0, (16) 

and the corresponding one for the other light quarks and their antiparticles. 
It allows to get the q(x) and Aq(x) distributions from the ones for q(x) and 
Aq(x). 

By performing a next-to-leading order QCD evolution of these parton 
distributions, we were able to obtain a good description of a large set of very 
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precise data on F%' n (x, Q 2 ), xF£ N (x, Q 2 ) and g{' d,n {x,Q 2 ), in correspon- 
dance with eight free parameters 15 . Therefore crucial tests will be provided 
by measuring flavor and spin asymmetries for antiquarks, for which we ex- 
pect 15 

Au{x) > > Ad(x) , (17) 

Au(x) - Ad(x) ~ d(x) - u(x) > . (18) 

For illustration, we show in Fig. 1 the predictions of the statistical approach 
with recent polarized DIS results. 




Figure 1. Results for A(u + u)/(u + u) (x) and A(d + d)/(d + d) (x) from Rcfs. [18,19], 
compared to the statistical model predictions. 

The inequality d{x) — u{x) > has the right sign to agree with the 
defect in the GSR 6 . Let us make a few comments on the BSR. In the low 
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x region ( x < 0.1 ), both g\ and g" are expected to increase, as well as 
their difference. The statistical prediction is close to the curve shown in 
Fig. 2, which has the simple expression 0.18(x -0 ' 5 — 1) implying 0.18 for the 
Bjorken integral. The statistical model gives the value 0.176, in excellent 
agreement with the QCD prediction 0.182 ±0.005 and with the world data 
0.176 ± 0.003 ± 0.07 20 . We also note that if Eq. (18) above is satisfied, 
it means that the antiquark polarization contributes to the BSR. In the 
statistical model this contribution is 0.022, which is not negligible. Finally 
this strong rise in the low x region, which was first noticed in Ref. [21] 
is consistent with the results from the resummation of double-logarithmic 
contributions 22 . 




Figure 2. The data on g\ n (x,Q 2 ) from Ref. [20] compared to the curve given by 
0.18(>- - 5 - 1). 
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4. The generalized GDH sum rule revisited 

The generalized GDH sum rules 14 are just being tested experimentally for 
proton, neutron and deuteron 23 ^ 24 ^ 5 ^ 6 . The characteristic feature of the 
proton data is the strong dependence on the four-momentum transfer Q 2 , 
for Q 2 < lGcV 2 , with a zero crossing for Q 2 ~ 200 - 250McV 2 , which 
is in complete agreement with our prediction 27 < 28 ; published almost ten 
years ago. Our approach is making use of the relation to the BCSR for 
the structure function g 2 , whose elastic contribution is the main source 
of a strong Q 2 -dependence, while the contribution of the other structure 
function, gx = gi + <?2 is smooth. However, the recently published proton 
JcffLab data 26 lie below the prediction, displaying quite a similar shape. 
Such a behaviour suggests, that the reason for the discrepancy may be the 
oversimplified treatment of the QCD expressions at the boundary point 
Qo ~ lGeV, defined in the smooth interpolation between large Q 2 and 
Q 2 = and which serve as an input for our model. For large Q 2 we took the 
asymptotic value for the GDH integral and we neglected all the calculable 
corrections, as well as the contribution of the g 2 structure function. This 
was quite natural and unnecessary ten years ago, since no data was available 
at that time. 

In a recent work 29 , we have filled up this gap by including the radiative 
(logarithmic) and power QCD corrections. The procedure used to take 
these corrections into account is explained in great details in Rcf . [29] . We 
found that the JeffLab data are quite sensitive to power corrections and may 
be used for the extraction of the relevant phcnomcnological parameters. In 
Fig. 3 we display for T\{Q 2 ), the comparison between the results of our 
previous work 27,28 and the new analysis which leads to a curve fairly close 
to the JeffLab data 26 . We show in Fig. 4 the same comparison for r™(Q 2 ) 
and we notice that the strong oscillation around Q 2 — lGcV 2 , we had 
in the previous analysis, is no longer there. These results can be easily 
transformed into predictions for the Bjorken integral T\~ n {Q 2 ) and for the 
deuteron Tf(Q 2 ) which turn out to be in good agreement with preliminary 
JeffLab data 30 ' 31 , in the low Q 2 region. 
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Figure 3. Our prediction for F^(Q 2 ). The thick line is the new analysis 29 , to be 
compared with the thin line, which represents our previous approach without corrections. 




Figure 4. Our prediction for r™(Q 2 ). The thick line is the new analysis 29 , to be 
compared with the thin line, which represents our previous approach without corrections. 
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